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Abstract 



Using the Operator Product Expansion for Wilson loops we derive a simple formula giving 
the discontinuities of the two loop result in terms of the one loop answer. We also argue 
that the knowledge of these discontinuities should be enough to fix the full two loop answer, 
for a general number of sides. We work this out explicitly for the case of the hexagon and 
rederive the known result. 
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1 Introduction 

Scattering amplitudes in N' = 4 super Yang Mills are related to polygonal Wilson loops with 
null edges [1]. The conformal symmetry on the Wilson loop side corresponds to non-trivial 
integrability charges on the amplitude side. For the case of ordinary operators, the conformal 
symmetry leads to the Operator Product Expansion. For Wilson loops we have a similar 
expansion, whose systematics was explained in [2]. 

In this paper we use this expansion to predict the form of certain discontinuities of the two 
loop expression for Wilson loop correlators. Namely, the Wilson loop correlators are func- 
tions which contain logarithmic branch cuts, corresponding to branch cuts at multi particle 
production thresholds for scattering amplitudes. Some of the discontinuities of the ampli- 
tude are determined by the operator product expansion. In the weak coupling perturbative 
expansion, this discontinuity can be obtained from lower order results plus the knowledge of 
certain anomalous dimensions. 

In the case of M^'^ kinematics we showed in [3] that we could determine the two loop results 
[4] from the one loop results plus the anomalous dimensions. In this paper we consider the 
general case of full M}'^ kinematics and we derive an expression for the discontinuities of the 
answer. We further conjecture that these are all the discontinuities of the two loop answer. 
The expression we derive for the discontinuity has the form of an integral kernel applied to 
the one loop result. We have not been able to do these integrals in general, however we have 
shown that one can rederive the known result [5, 6, 7] for the hexagon using this method. 

The discontinuities are particularly simple at two loops. The reason is that the one loop 
result comes from the tree level exchange of a gluon between different sides of a Wilson loop. 
Thus, the OPE at one loop contains only single particle exchanges. These particles behave 
as free particles and transform under simple representations of the full conformal group. At 
two loops we get a logarithmic term in the OPE which comes from the anomalous dimension 
of these states. These anomalous dimensions were computed in [8]. This logarithmic term 
gives rise to a discontinuity of the two loop answer. In fact, if we include all states in 
the OPE expansion, we can multiply each term by the corresponding anomalous dimension 
and perform the whole sum. This gives us an expression for the discontinuity. This whole 
operation reduces to a simple kernel acting on the one loop answer. This is one of the main 
results of the paper. 

In addition, we observe that these discontinuities appear to be enough to determine the 
answer. This can be most simply explained by using the "symbol" operation introduced 
in [7]. Namely the functions appearing in the two loop answer can be assigned a symbol 
which is a sum of terms of the form Sym[yl] = ^ -Ri ® i?2 ® -R3 ® -^4- We conjecture 
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that the first term in the symbol is always a distance between two cusps. This is simply 
saying that we expect a discontinuity whenever two cusps, which were originally spacelike 
separated, become null separated. The number of independent discontinuities is the same 
as the number of independent OPE expansion channels that we can have. Our expressions 
for the discontinuity are transcendentality three functions which determine the last three 
factors of the symbol. Thus, if we had an efficient way to determine the symbol of our 
integral expression, we could determine the symbol of the answer. 

As we mentioned above, we have carried out this program for the simplest case of the 
hexagon. In this way we have rederived the known answer for the hexagon. Of course, in 
doing so, we have explicitly shown that the known two loop expression has a proper OPE 
expansion, in agreement with the general discussion in [2]. 

This paper is organized as follows. 

In section two, we give a brief review of symbols and their properties since we found them 
useful for stating our results and also for deriving some of them. We also explain why the 
various OPE channels determine the symbol of the two loop answer, if one assumes that the 
first entry in the symbol is a distance. In section three we compute the OPE expansion at 
one loop. Namely, we write the part of the one loop answer that contributes to the OPE. 
This is given by a correlation function of two Wilson loops. We then decompose this answer 
into the contribution from self dual and anti-self dual fields. This is necessary since the 
anomalous dimensions will treat these two cases differently. 

In section four we describe the various states that are propagating in the one loop result. 
We write their anomalous dimensions, originally derived in [8]. We also show that we can 
describe the action of the anomalous dimension as certain convolution kernel. 

In section five we describe the hexagon. We perform the expansions explicitly and we 
compute the action of the convolution kernel indirectly, by considering the symbol and the 
possible discontinuities. This is a method that could be useful in the general case, though 
here we have only derived it for the case of the hexagon. 

Some technical details are relegated to an appendix. In addition, we have two appendices 
with two special kinematic limits where the answer simplifies and can be computed more 
easily. 

2 The OPE and the Symbols 

We have found that the notion of a "symbol" introduced in [7] (and references therein) is 
very useful for dealing with the type of transcendental functions that appear in these Wilson 
loop computations. Our final answer can be stated without any reference to the symbol. 
However, it was useful to use symbols for intermediate steps, as well as for thinking about 
these functions. Thus we will begin with a quick review of symbols. 
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2.1 Introduction to symbols 



Transcendental functions such as polylogarithms and their generahzations satisfy many in- 
tricate identities. This imphes that given two combinations of polylogarithms, it can be very 
tricky to recognize if they coincide or not. There is a linear operation on transcendental 
functions, which we will denote as a map F — )■ Sym[F], which is very useful to keep track 
of all sort of polylogarithm identities. The polylogarithm identities satisfied by a function 
F become trivial algebraic identities satisfied by its "symbol" Sym[F]. The symbol Sym[F] 
loses some of the information about F, but it is still very useful. 

Functions of transcendentality degree n are defined recursively, as functions with loga- 
rithmic cuts, such that the discontinuity across the cuts is 27ri times functions of transcen- 
dentality degree n — 1. Constants (or rational functions), by definition, have degree 0. We 

can give some simple examples. The logarithm log a; has degree 1: is has a cut running along 
the negative real axis, with discontinuity 27ri. The dilogarithm Li2{x) has degree 2: It has a 
cut running from 1 to oo with discontinuity —27[i\ogx. 

The symbol of a function of transcendentality n is a linear combination of elements of 
the form 

Ri ^ R2 ^ Rs ^ Ra - ■ ■ ^ Rn (1) 
where the Ri are rational functions. It should be thought of as a tensor product of logarithms 

log Ri ^ log R2 log i?3 log -R4 • • • ® log Rn (2) 

but the "log" is usually omitted. We have the following recursive definition: 

Sym[F] = J^i?„®Sym[Fj (3) 

a 

if F has logarithmic cuts which start at Ra = 0, end at Ra = 00, with discontinuity 27iiFa. 

For example, Sym[l] = 0, Sym[loga;] = x, S'ym[—Li2{x)] = {x — 1) ^ x. The symbol 
of log a; log 1/ is x <S) y + y <S) x: it has discontinuity 27rilog|/ around x — and 27rilogx 
around y — 0. Notice that the difference of two functions of degree n with the same symbol 
can only be a function which becomes after taking n discontinuities, i.e. a function of 
transcendentality degree n — 1 or lower. This is the amount of information which is lost by 
taking the symbol. 

From the definition, a lot of useful properties follow right away. The most important are 
that the symbol behaves like a tensor product of logarithms 

■■■®xy0--- ^■■■®x® \ (g)y(8)--- (4) 

and it is transparent to constants: 

• • • ® (ci?i) ® ■■■ = ■■■ ® i?i ® •• • (5) 

if c is a constant. We can illustrate this for a simple dilogarithm identity: 

TT^ 1 

Li2{x) + Li2{l/x) = -— - - log(-x)2 
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becomes, upon taking the symbol, 



(1 - x) (g) X - (1 - 1/x) (g) (1/x) 



(1 — x) (g X + 



X — 1 



(g X = ( — 1/x) ® X = —X ® X 



X 



The symbol of the product of two functions Sym[FG] can be built from the relation on 
discontinuities A{FG) = {AF)G + FAG. By linearity of the tensor product, we only need 
to compute it for Sym[F] = (g)"=i-Ri and Sym[G] = (g^^+i-Rj: 



where the sum is over all the permutations of the full set oi n + m Ri which preserve the 
original ordering among the factors in Sym[F] and among the factors in Sym[G]. 

For example, if n = m = 2, 

Sym[FG] =Ri ® R2 ® R3 ® R4 + Ri ® R3 ® R2 ® Ra + Ri ® R3 ® Ra® R2+ (7) 
+ R3(^ Ri® R2® R4 + R3® Ri® Ra® R2 + R3® Ra® Ri® R2 (8) 

Hence we can easily build the symbol of any polynomial in polylogarithms and logarithms. 

There is a relation between the symbol of a function and the definition through iterated 
integrals. Indeed, the operation of taking a discontinuity commutes with differentiation, 
hence the symbol of a derivative Sym[(iF] is given by the replacement in Sym[F] 



Not all possible symbols are symbols of a function. For example, if x and y are unrelated, 
x^y is not the symbol of a function. A useful test to check if a symbol could be the symbol 
of a function is to pick two consecutive slots of the symbol, and replace [9] 

-Ri ® ■ ■ ■ ® -Ri (g Ri+i (g ■ • ■ ® -R„ -> dlogRi A (ilogi?i+i ® ■ ■ ■ » • • • ® i?„] (10) 

This transformation maps the symbol of functions to zero. 

2.2 On OPE constructibility 

The null polygonal Wilson loops which we study are defined as the limit of Wilson loops with 
space-like sides, and space-like separation between any pair of points. In the limit, we make 
the sides of the polygon null, but focus on the kinematic region where all points which are 
not on the same side remain space-like separated. Infrared divergencies which arise in the 
null limit are well understood. The answer should be smooth in this kinematic region. The 
kinematic region of interest is bounded by walls where two non-consecutive vertices become 
null-separated, i.e. (xj — XjY = 0. We can analytically continue the amplitude around these 
loci, but the analytically continued answer will generically have cuts which end at these loci. 




(6) 



TT 



i?i ® ■ ■ ■ ® i?„ dlogRn [Ri®---® Rn^l] 



(9) 
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At one loop, the Wilson loop amplitude is a function of transcendentality degree two. It 
has logarithmic cuts which start at (xj — xj)'^ = 0. More precisely, if we take the symbol of 
the one-loop answer, it basically takes the form 

51-ioop ^ ^(^^ _ ^ Sjm[Dl7'^°^] (11) 

where dI~^"°^ are the discontinuities at the corresponding cuts. 

The two-loop answer will in general also have cuts which start at the loci (xj — Xj)"^ = 0. 
In order to bootstrap the two-loop amplitude, we will make a crucial assumption: we will 
require the two loop answer to be a function of transcendentality degree four, which has 
transcendentality degree three discontinuities only across the cuts starting at the the loci 
{xi — XjY = 0. More precisely, we will require the symbol of the two- loop answer to take the 
form 

^2-ioop ^ ^(^^ _ ^ Sym[Dj,-^°°P] (12) 

We will now show how to compute the discontinuities around the loci (xj — XjY = by 
an analysis of OPEs and hence compute the full symbol of the two loop amplitude. Notice 
that a generic formal symbol is not the symbol of an actual function. Given the set of 
Sym[Dj?~^°°^], we can check if 5*^"'°°^ in (12) is the symbol of a function or not. If it is, 
our basic assumption passes a very strong self-consistency check. If it is not, it means that 
the two loops amplitude does not have a symbol or alternatively, that our original set of 
discontinuities was incomplete, and the amplitude has other cuts further away from the 
original kinematic region, which need to be accounted for in the boot-strap procedure. 

2.2.1 Discontinuities from OPEs 

The main idea behind OPEs is to study (multi)collinear limits of the amplitude, by acting 
on a portion of the Wilson loop with a one-parameter family of conformal transformations 
Mr, chosen in such a way that as r — )■ 00 that portion of the Wilson loop is pushed towards 
a (multi) coUinear limit [2] . The distances between the vertices in that portion of the Wilson 
loop scale uniformly as for large r. The family of Wilson loops remains inside the 
preferred Euclidean kinematic region, and approaches the boundary as r — )■ cxd. An OPE 
expansion channel is chosen by selecting a pair of null sides separated by two or more other 
null sides. We call one of the sides the "left" side and the other the "right" side. 

The general analysis in [2] constrains the behavior of the amplitude at large r by an OPE 
analysis, which replaces the portion of Wilson loop by a sum of operator insertions on the 
T ^ 00 limit of the Wilson loop. At two loops, the result is an expansion for the amplitude 
of the form 

60 + {aiT + bi)e-^ + (aaT + 62)6"^^ + (agr + 63)6"^^ + ■ ■ ■ (13) 

If we complexify r, we can choose a path in the space of Wilson loops in order to compute 
useful discontinuities: go to very large r, then take the discontinuity of the amplitude under 
the shift r — > r + 27ri, which is a small closed loop in the space of Wilson loops. Upon 
taking this discontinuity, only the re"""^ terms contribute. These are very special: they arise 
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from the "anomalous dimension" of the operators, i.e. the leading correction to their quan- 
tum numbers under the Mr transformation. Hence the OPE discontinuity of the two-loop 
amplitude differs from the full one-loop amplitude only by the insertion of this anomalous 
dimension operator. 

This path at large r, winds once around all the loci (xj — XjY = for all pairs of vertices 
i, j in the portion of the path acted upon by Mr- This "OPE discontinuity" is 

where the sum runs over all pairs of vertices z, j in the portion of the path acted upon by 
Mt-. These are the vertices which are between the left and right sides as defined above. 

This means that we can reconstruct all the recursively. At the first step, we 

extract -0^7+2°'' from the OPE discontinuity corresponding to the collinear limit of the 
i,i + l^i + 2 vertices. At the second step, we extract -0^^7+3°^ from the OPE discontinu- 
ity corresponding to the collinear limit of the i,2-|-l,i-|-2,i-|-3 vertices, which is 

Dl,r = + itS" + D^-'^l, (15) 

and so on. 

Thus, if we have a way to computing all the OPE discontinuities for all OPE channels, 
namely, for any choice of two sides L, i?, we can then construct the symbol for the two loop 
answer (12). Below we will derive an expression for this discontinuity, as a function. In 
fact, Dl,r, and thus, also D^~^°°^ are symbols of a function. This also follows from the fact 
that we can vary the Xi.Xj independently. Note that conformal invariance implies relations 
among the Df^^"°^ = Z}^7^°°^) ^ind J2i d1j^°"^ = for each j. 



3 The one loop sources 

In order to compute the two loop discontinuities we need to compute the one loop contri- 
bution to the OPE. This is simply the one loop result for the Wilson loop where we keep 
only the part that depends non-trivially on the OPE "time" parameter r. This is completely 
captured by a single particle exchange between the "top " part of the polygon and the "bot- 
tom" part. The top and bottom parts are separated by the two special sides, called "left" 
and "right". Selecting and OPE channel is the same as selecting these two special sides. This 
can be derived from [10]. We rederive these results from scratch here. 



3.1 Short review on kinematics 

We represent M}'^ as the lightcone in M^'^. That is, any point x G M}'^ is associated with a 
null ray in M^'^, {X G M^'^, = 0, X ~ tX}. The map to the usual Poincare coordinates 
is 
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The conformal group is then reahzed hnearly as the SO (2, 4) rotations of the embedding 
coordinates (the X's). A convenient way of representing a null polygon is using momentum 
twistors [11]. These are M^'"^ spinors that are defined up to rescaling A = (Ai, A2, A3, A4) ~ tX. 
The conformal group acts on these by a simple >S'L(4) left group multiplication. With any 
pair of momentum twistors A, A we can associate a null ray in M^'^ and therefore a point 



X e M^'^ as 



V —V _ \ \ M — — 1, 0, 1, 2, 3, 4 , , 

y^ab — Am! ab — ^[a^6] , a b = 1 2 3 4 ' ^ ^ 

where Tm are sigma matrices in M^'^. A polygon with null sides can be given as a sequence 
of A^ twistors A,, such that the intersection of the sides i and i + 1 is the point Xj = Aj A Aj+i. 
The distance between two cusps of the polygon reads 



(Xi - X 



_ Aj A Aj_|_i A Xj A Aj+i 

" (Af ■r+-A,+i)(Aj-r+-A,+o 



where T"*" = (F ^ + F^). The numerator in (18) denotes the SL(4) invariant product of 
four twistors, given by contracting with the e°'^"'°"^ tensor of SL{4). We automatically have 
Xi_i ■ Xj = 0, which is the condition that the i-th side should be null. It is also useful 
to introduce dual momentum twistors Aj, which can be defined as Aj = Ai_i A Aj A Aj+i. 
They also satisfy Xj = Aj A Aj+i and transform by a right 5*^(4) multiplication. A twistor 
A together with an orthogonal dual twistor fi, {fi, A) = Aa/t" = 0, defines a null line (and 
viceversa). The line is given by all points of the form Z = X Av with (/t, v) = 0. 



3.2 One-loop pairing 

In this section we compute the integral of the photon propagator stretched between a pair 
of closed null-polygonal Wilson loops with relative space-like separations (see figure 1). 

The calculation is in principle straightforward. We will use a six-dimensional description 
of conformal space-time. Points in spacetime correspond to equivalence classes of null six- 
dimensional vectors X up to rescaling. A possible choice of representative takes the form 
X = (1, x^, x"^), so that X -Y = {x — y)"^. We will pick two six dimensional representatives 
X(s), Y{s') for the paths. Then the integral 

f fdX-dV 

r^(i)(X,r) = y y (19) 

will do the job. The rescaling invariance of the integral is a bit subtle: a rescaling X(s) — ?■ 
A(s)X(s) will shift the integrand by dsXdg' logX ■ Y. As long as the polygons are space-like 
separated, so that the logX ■ Y has no monodromy, the integral is conformal invariant. In 
the conformal gauge where X ■ Y coincides with the four-dimensional squared distance, the 
integrand "^^'^^ reduces to the desired j^i^- 

The integral can be done by brute force, say by picking the parameterization X(s) = 
(1 — s)Xj + sXj_|_i in the i-th side of the first polygon, and similarly for the second polygon. 
The result of the side-by-side brute force calculation is a sum of dilogarithms and bilinears of 
logarithms, but the arguments of the dilogarithms contain linear combinations of the Xj ■ Yj 
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(a) (b) 

Figure 1: Correlation function between two Polygonal Wilson loops, which is given by a 
single gluon exchange. In (a) we consider a Polygon where all the points in the X contour 
are spacelike separated from all points in the Y contour, (b) We take a limit of (a) where a 
segment of one contour {Y^, Ife+i) is on the same null hne as a cusp of the other contour Xj. 
The rest of the points are spacelike separated. 



which are not well-behaved under rescaling. Dilogarithm identities need to be used in order 
to eliminate this spurious dependence. A shortcut to deal with such identities is to compute 
the symbol of the answer. We refer to section 2.1 for a review of symbols of poly logarithms. 

The integral {X, Y) has interesting discontinuities around the loci where the polygons 
become null separated first. It is not difficult to see that generically this happens first 
when a vertex of a polygon becomes null separated from a vertex in the other polygon. 
Indeed if a smooth point Xq = X{sq) for X{s) is null separated from a point Yq of Y{s'), 
X{s).(Yo + 6Y) ~ (s — so)X(so) ■ Yq + Xq ■ 6Y so either a whole side of X is becoming null 
separated to Yq or we can find a solution for s — sq describing points of X near Xq null 
separated from points of Y near Iq- The symbol of the full answer turns out to have the 
form 

5- = Sym[r^(i)(X,r)] = J^X, ■ ® Sym[A,fc]. (20) 



i.k 



loci, with a 



The integral has a logarithmic cut originating from each of the Xi ■ Y^ 
discontinuity 27riAjfc. 

If two vertices Xi, Yk are close to be null separated, one can study directly the local 
behavior of the integral near the two vertices. The result is that 



Aifc = log 
-log 



1 - 



(x,-n)(x,_i-n+i) 

(x,-n+i)(x,_i-n) 

i^i ■ ^/c-l)(-^i-l ■ ^k) 



+ log 
- log 



1 - 



(X,-n)(X,+i-F, 



k-l 



(x,-n_i)(x,+i-n) 

i^i ■ ^fc+l)(-^j+l ■ ^k) 



(21) 



10 



We can now write down a function with symbol S, and no discontinuities in the region 
where the polygons have space-like separation. This should agree with r[/(i) (X, Y) up to a 
constant, which can be fixed simply by requiring the answer to go to zero as the two polygons 
are brought far from each other. 

We can rewrite S as 



5 = E 



ik 



i^i ■ yk+l){^i-l ■ ^k] 



i^i ■ ^fc+l)(-^i-l ■ ^k] 



■ Yfc_i)(Xj_i ■ Yk) (Xj ■ Ffc+i)(Xj_|_i ■ Yk] 
The first line is the symbol of a sum of dilogarithms 

■ Yk){Xi_i ■ Yk+i] 

ik 



ri 



(22) 



(23) 



and the second part can be reorganized a bit to the symbol of the rescaling- invariant com- 
bination 

E, V J (Xi ■ Yk){Xi+i ■ Yk+i) 

Indeed a rescaling of Xj or of Yk shifts r2 by a sum which telescopes to zero. Then 
r(/(i)(X, y) = ri + r2 in (23) and (24). Notice that in the kinematical region of space- 
like separation, we can take Xj ■ > 0, and the functions above have manifestly no cuts. 
Notice also that the discontinuity function Aj^ are zero at the boundary Xj ■ = 0, which 
means that the ^^/(i) (X, Y) does not diverge if a vertex of X and a vertex of Y are brought 
to be null-separated. 

There is a useful rearrangement of the answer in terms of momentum twistors. Mo- 
mentum twistors are four-components vectors which transform in the fundamental of the 
SL{4) conformal group. A six-dimensional null- vector X can always be written as the 
wedge product of two momentum twistors (17). For a generic pair of polygons, we can 
pick momentum twistors Aj associated to sides of X and /i^ associated to sides of y, so that 
Xj = Aj A Aj+i and Yk = Hk /\ jJ^k+i- Then inner products are written as determinants, 
Xi ■ Yk = (Aj A Aj+i A /ifc A /ifc+i) and we have a very useful Plucker identity: starting for 
example with 

{Xi ■ Yfe+i)(Xj_i ■ Yfc) — (Xj ■ Yfc)(Xj_i ■ Yk+i) = (25) 

(Aj A Aj+i A fik+i A /ifc+2)(Aj-i A Aj A yUfc A yUfc+i) 
-(Aj A Aj+i A /ifc A Aifc+i)(Aj_i A Aj A yUfe+i A /ife+2) 

we observe that Aj A Afc+i is common to all determinants. Modulo Aj, A^+i we can use a 
standard two-dimensional cyclic identity and rewrite 

{Xi ■ Yfc+i)(Xj„i • Yfc) — (Xj ■ Ffc)(Xj_i ■ Yk+i) = (Aj A /i^ A fik+i A fik+2){f^k+i A Aj_i A Aj A Aj+i) 

= (Ai,/ifc+i)(yUfc+i, Aj) (26) 
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where we have introduced the dual momentum twistors for the contours, as Aj = Aj_i A A,, A 
Aj+i and fik+i = fik A fik+i A fik+2- We can combine the four logs in Aj^, the denominators 
drop out and we get simply the product of two twistor cross-ratios 



log 



(Ai, /ifc+i)(Ai+i, /ifc) (/ife+i, Aj)(/ifc, Aj+i) 
(Aj, /ife)(Ai+i, fck+i) (/ifc, Aj)(yUfc+i, Ai+i) 



(27) 



This is a very neat answer. Notice that the symbol 5* of ?";7(i) {X, Y) can be now split into 
two pieces: 

(Ai, /ifc+i)(Aj+i, fik) 



and 



s- 



(Ai, /ifc)(Ai+l, /ifc+l) 

ifJ'k+l, ^i){fJ'k, Ai+l) 
(/^fcj Ai)(/ifc+i, Ai+i) 



Xi ■ n ® A, 



ifc 



Xi ■ Ffc ® A, 



ik 



{2i 



(29) 



We can even write down separate functions with these two symbol. The crucial 
observation is that 

Xi-Yk = (Aj, fik){K+i, fi'k+i) — (Ai, fi-k+i i)Afc) — (/^fc? Ai)(/ifc+i, Ai+i) — (/ifc+i, Ai)(/ifc, Aj+i) 
Hence we can subtract off from 



i,k 



1 - 



(Ai, /ifc+i)(Ai+i, fik) 



(Ai, /ifc+i)(Ai+i, /ifc) 



(Ai, /ifc)(Ai+i, ^k+i) _ 

which is the symbol of 

^ .(Aj, /iA;)(Ai+i, /ifc+ij_ 

and be left with the symbol of a bunch of logarithms 



(Ai, /ifc)(Ai+i, /ifc+i) 

(Aj, flk+\){\+li P'k) 



^log(Ai,/ifc)log 



i.k 



(Ai, /ifc+i)(Ai+i, /tfc) 
(Ai, /ifc)(Ai+i, /ifc+i) 



(30) 



(31) 



(32) 



Thus = + in (31), (32). r~ is given by the same expressions exchanging twistors 
and co-twistors, namely Ai — )■ Aj, fik — )• fJ.k- And finally, r[/(i) = r+ + r~. Actually, if we 
want to have answers which explicitly have no cuts in the region of space-like separation of 
the two loops we should pick the arguments of the dilogarithms with a little more care. The 
boundaries Xi ■ = are the locus where the argument of the corresponding dilogarithm 
becomes 1. In a region of space-like separation, some arguments will be strictly larger than 
1 and some strictly smaller than 1. Indeed the product of the arguments over k at fixed i 
telescopes to 1. In any component of the region of space- like separation we should do some 
replacements Lz2(l/x) — ?■ —Li2{x) — ^log^a; so that all dilogarithms have no cuts on the 
correct side of 1. 

For our computation of the two loops remainder function we have to consider the case 
where a cusp of the bottom polygon Xi is null separated from an edge of the top polygon 
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(FfcjYfc+i) and a cusp of the top polygon Yj is null separated from an edge of the bottom 
polygon {Xi,Xi+i), see figure 1(b). In this case, reduce to manifestly finite functions. 
This can be found explicitly in appendix C. 

In order to understand better the meaning of and the corresponding functions r^{X, Y), 
it is useful to review a bit of twistorial geometry. 

3.3 Null planes and null lines 

In (2, 2) signature it is possible to find closed null triangles. Given three momentum twistors 
•^i) -^2 5 -^3 we can make a triangle the standard way, with vertices Xj = Aj A Xi+i- The whole 
triangle satisfies the constraint X ■ X = 0, with the dual momentum twistor A defined as 
A = Ai A A2 A A3. This is a rather interesting equation. It defines a null plane, which we 
will denote as 5*. Two such null planes either coincide, or intersect at a single point: the 
equations A ■ X = and A' ■ X = have a unique solution (up to rescaling) for X. 

There is a second type of null planes, where we exchange the roles of twistors and dual 
twistors. They are planes containing the second type of possible null triangles: triangles 
with vertices of the form Xj = A A Aj for some A. The points X = A A f for some fixed A, 
variable v generate a null plane which we will indicate as S. We can also describe it as the 
space of solutions to the equation X A A = 0. Again, two null planes S and S' either coincide 
or meet at a unique point, X = A A A'. 

A S plane and a S plane with A ■ A 7^ will not intersect at all: A ■ (A A f ) is never zero. 
On the other hand, if A ■ A = the two planes meet along a null line: the set of points such 
that X = X Av with X ■ v = 0. Vice- versa, any null line can be described this way, in terms 
of a twistor and a dual twistor: given two null-separated points Xi and X2 we can always 
write Xi = A A Ai and X2 = A A A2 and then define a A = A A Ai A A2. The two points and 
any linear combinations of them lie on the line defined by A, A. So a pair of null-separated 
points define two null planes of opposite type and a null line through them. 

Given a null plane S and a point X not on it, the twistor A = A • X defines a second 
null plane S which passes through X and intersects S along a null line. The same is true for 
a null plane S and a point X not on it, though then it is A = X A A that defines a second 
null plane S and null line. Given a line i parameterized by a pair (A, A) and a point X 
away from the line, we can find a unique point on the line which is null separated from X: 
X^ = (A ■ X) A A. The null line through X and Xi is given by the pair (A ■ X, X A A). 

3.4 Helicity interpretation of r^(X, y) 

It is interesting to compute the pairing ^[/(i) (X, Y) of two null triangles of the same type. If 
we look at two triangles of the same type, say with vertices Xj = A A Aj and = A A A^, 
we notice that any inner product Xj ■ 1^ = A A Aj A A A A^ involves the quantity AAA, 
hence depends on Xi,Xk only through their value modulo A, A. This reduces the cross-ratios 
effectively to 2d cross-ratios, and it is easy to check that Ajfc = 0. Hence r(7(i)(X, F) (or 
at least its symbol) vanishes for two triangles of the same type. This has a simple physical 
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Ai A A2 A2 A A3 Ai A A2 A2 A A3 



Figure 2: Triangulation of a null polygon Wilson loops into the two type of triangles. The 
two types of triangles correspond to the ones bounded by A's which are all different or the 
ones with the same A for all three edges. 

interpretation: the integral of the U{1) gauge field along a triangular contour is the same 
as the integral of the field strength on a surface bounded by the triangle. If the triangle is 
null, it lies on a null plane and we can take the surface to lie on the plane as well. But the 
surface elements of the two types of null-planes are self-dual and anti-self-dual respectively, 
and hence couple to photons of positive or negative helicity only. But the propagator is 
non-zero only for zero total helicity. 

The pairing ^[/(i) (X, Y) for triangles of opposite type is a non-trivial function of the 
three independent cross ratios which can be built from the nine inner products Aj ■ /t^ of the 
three twistors which define a triangle, and the three dual twistors which define the other 
(or vice- versa, depending on the type of triangles). The symbol S can be computed as done 
for general polygons, but it differs on one important way: it only includes one of the two 
twistorial cross- ratios, hence S = S"^ {or S = S~ , depending on the type of triangles). 

This is possibly interesting, because any null polygon can be decomposed into a sum of 
triangles of the two types. Consider any three consecutive sides of a polygon, associated with 
twistors Ai, A2, A3, and imagine replacing the sequence of three sides with two new sides, 
defined by the new vertex Ai A A3. It is easy to see that this point is null-separated from the 
initial and final vertices of the sequence of sides. The new polygon has the same sequence 
of twistors as the old one, with A2 omitted (see figure 2). The difference between the old 
polygon and the new polygon is a pentagon. But the new vertex is actually null-separated 
from all points in the old sequence of sides. The pentagon is really the sum of three null 
triangles, one of the first type (lying in the null plane Ai A A2 A A3) and two of the second type 
(lying in the null planes Ai and A3). Of course, it is possible to invert the roles of twistors 
and dual twistors in the above operation. 

Iterating either operation, in any order, one can produce a variety of triangulations for 
any given polygon into triangles whose vertices are a subset of all pairs of twistors, or all 
pairs of dual twistors in the polygon. Geometrically, these are various intersections of the 
null-planes passing through the sides of the polygon. 

Now, pick any such triangulation of two space-like separated polygons and write r(/(i) (X, Y) 
as a sum over pairs of triangles in the two triangulations. We can naturally decompose 
ru{i) {X, Y) in two pieces: r"*" includes the contributions from triangles of the first type in X 
paired up with triangles of the second type inY, r~ includes the contributions from triangles 
of the second type in X paired up with triangles of the first type in Y. 
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Notice that the symbols of contain only either type of twistorial cross-ratio, but have 
to add up to the symbol S of ru(i){X^Y). Hence they must coincide with and S~ , and 
do not depend on the choice of triangulation. The actual functions do not depend on 
the triangulation either: we can simply focus on one triangle in X, say of the first type, and 
sum over the triangles in Y . This only contributes to r"*", but it is the full r{7(i) between the 
triangle and Y , hence it is independent of the triangulation of Y . Summing over triangles of 
the first type in X we confirm that r"*" is independent of the triangulation of Y . The same 
reasoning with triangles of the second type in Y proves independence on the triangulation 
oiX. 



3.5 Differential equations for r (X, F) 

In the OPE analysis, we will need to decompose ru{i){X,Y) as a sum of contributions 
corresponding to the propagation of a particle with given quantum numbers under various 
choices of conformal transformations. In other words, we will need to decompose ru(^i)(g o 
X, Y) in pieces which depend in specific ways on the SL^A) conformal transformation matrix 
g. We will particularly interested in the decomposition under a certain SL{2)^ x SL{2)r x 
U{1)^ subgroup of SL{4:) 

fg^e'^'^ \ 
^=^0 g.e-^^l-) (33) 

These transformations preserve an M}'^ inside E?'^: U (1)^ rotations around the M}'^^ SL(2)cr x 
SL{2)j. conformal transformations in the M}'^ plane. Alternatively, they are the isometrics 
of a AdS^ ® spacetime, conformally equivalent to flat space. 

When acting with g on the X oi Y polygons, we find that the two helicities r+ and r~ 
give rise to two towers of excitations, with a specific relation between the r and a spins and 
the U{l),j, quantum number m. We can check directly the consequences on ru{i){g o X,Y) 
of this fact. If we take Aj to be in the fundamental representation of g, the Casimirs Co- and 
Cr of SL{2)^ and SL{2)r satisfy 

when acting on r~^{X, Y), and 

a = ^ + y (35) 

when acting on r~(X, Y). These Casimirs are given by s{s — 1) and we can call s the spin. 
Note that s and 1 — s give the same Casimir. 



4 Convolution with the anomalous dimension kernel 
4.1 The anomalous dimension 

In this section we recall the formula for the anomalous dimension of the operators that we 
are exchanging. 
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As we discussed, in our problem it is convenient to focus on an SL{2)^ x SL{2)„ x S0{2) 
subgroup of 5'L(4) [2]. We consider a U{1) gauge field and expand it into representations of 
these symmetries. The states are multiplets of SL{2)r and all states in a given multiplet have 
the same anomalous dimension [3]. The dimension is the Lq generator inside SL[2)r. The 
anomalous dimension depends on the momentum of the field. The momentum is a generator 
inside SL{2)„. It turns out that the anomalous dimension depends also on the spin under 
SL{2), of the field. 

Thus, it is convenient to classify the states into primaries under SL{2)j. and SL{2)„. It 
is convenient to think of the states as operators acting on the origin. In that language, the 
L_i generators of SL{2)r is d- and the corresponding one in SL{2)a is d+. The other two 
spacetime derivatives are just dz,dz, and carry charges ±1 under 50(2). The states are 
constructed from the field strength F^^, and its derivatives. For the purposes of classifying 
the states we can view the field strength as abelian and the derivatives as ordinary partial 
derivatives. We also should impose the equations of motion and Bianchi identities, dF = 
d * F = 0. We are then interested in states which cannot be written as d± of other states. 
We can consider the twist operator which is the dimension minus the spin in the H — plane, 
r = A — M_|__, where M+_ is the spin in the H — plane. We can define the conformal spin 
under SL{2)t as /3 = |. This assigns zero twist to The twist corresponds to the dr 
generator in SL{2)r- The states with lowest twist are the twist one states -F-g which 
have S0{2) charges plus or minus one. In this fashion we can classify all states and we get 
the following ones 



F+- 


— F - 


^z ^ z- , 


m > 1 




- 1 m < —1 


(36) 


F+- 


- + Fz, , 




m > 1 


d:^-'Fz^ 


- 1 m < —1 


(37) 




/3 = 1 , 






s = /3- 




(38) 


s = 


s- 1 = 


^ 2 ' 





We have indicated the conformal spins under SL{2)fj and SL{2)r in the last line (38). We 
see that we have two states at twist two, given by the first column. Then for each SO (2) 
charge m, m 7^ 0, we have a states with twists \m\ or \m\ + 2 given by the second and third 
column respectively. The conformal spin s of each state is given by s = ^^^^+- . 

The anomalous dimension of these states is then [8] 

72s(p) = 2/ [i,{s + i^-) + iPis - z|) - 2^(1)] (39) 

where if) is Euler's psi function ■ip{x) = Note that the first hne (36) involves self dual 

fields and the second line (37) involves anti-self dual fields. The conformal spin of the two 
towers is related to the S0{2) charge. Thus, we can rewrite (39) in terms of the S0{2) 
charges as 

For (36): 7-m(p) = 72+™(p) (40) 

For (37) : 7„(p) = 72-m(p) (41) 

Thus, we have now expressed the formulas for the anomalous dimensions in terms of the 
eigenvalues under some of the generators of the conformal group. Here m denotes the 5*0(2) 
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charge, including its sign. Note that if s is a half integer, then (39) is invariant under 
s — 7- 1 — s, which implies that the anomalous dimensions depend only on the casimir of 
SL{2)rj (and the momentum, which is an element in SL{2)cr)- 



4.2 The anomalous dimension as a convolution 

If we start from the one loop sources computed in the previous section and we Fourier 
transform them we can multiply each term by the corresponding anomalous dimension. 

Alternatively, we can transform a product of fourier transforms into a convolution. For 
that purpose it is useful to use the integral expression for the ip function 

oo 

^(1 + z)- ^{1) = -2jdt —^r^ (42) 



This integral expression is good for z > —1. Thus, in writing (40) we have to separate the 
cases m > and m < 0. We will get insertions of e^^^i^^i^^J) ^ where is the momentum 
in the sigma direction and J is the SO (2) charge, whose eigenvalue is m. When we insert 
this in the one loop source, it is the same as deforming the upper part of the polygon by an 
element of the conformal group. Using the integral expression we get that the convolution 
of the kernel with r+ ends up giving us an expression of the form 

7 rr+(e*(^^-^)X, Y) + r+(e*(-^^--^)X, Y) - 2r+(X, Y)] 
7*r+ = -Ag^jdt'- ^^^^ ^ (43) 







, 7 re2*{r+(e*(^^+-^)X,r) +r+(e*(-^^+^)X,r)| -2r+(X,r,, , 

I J J- '- ^ ' \ U 



-Ag^ / dt 



o2t 







Where the symbols [■■■]>o and [■■■]<o indicate that we project on to the part that has 
positive (or zero) eigenvalues (or frequencies) under J or negative eigenvalues respectively. 
In performing these projections we assume that we are near the OPE limit, which gives a 
natural size to various terms, or a choice of contour. The last term, the one that is not 
shifted by t can be dropped if one is computing the symbol, since it simply gives the one 
loop answer times a constant. 

These projection conditions can be done by performing contour integrals 

r>o{,z) = I r{w) , r<o{z) = --^ / -^^^r{w) (44) 

Zm J [w — zj 2m J w — z 

\w\>\z\ \M<\z\ 

When we use these expressions in the integrand of (43) we can view e^*"^ as acting on z. In 
other words, let us say that the original value of z, zo of the original polygon can be viewed 
as a point on a circle, Zq = e*"^. Then in the positive frequency projector (44), which acts on 
the first line of (43), we set z = e~*e*'^ and w = e*'^"'"*''^, and and integrate over (p. On the 
negative frequency projector we simply have z = e^e^'^ and the same w. In the end, we see 
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that r_|_ is evaluated at 2;oe"^. This is the same as saying that the X polygon is acted upon 
by e*"^"^ where J is the 5*0(2) generator. 

We can combine the two integrals and write the full kernel acting on the positive helicity 
contribution as 

oo 27r 

' ^ ^ J J 27r (1 - e-*e-*^)(l - e*e-*^) ^ ^ 

-oo 

The kernel for 7 * r_ is the same except that we replace e"*"^ — )■ e**^ in the kernel in (45). 
Alternatively, we can use an expression for 7 * r_ which is similar to (43) except that we 
exchange J -J and [■]>o [■]<o, [■]<o ^ [■]>o- 

It is possible to do one of the integrals in a straightforward way (using mathematica, 
for example). The second integral then looks a bit ugly. (It can be done but it involves 
some square roots). Perhaps it becomes easier after simplifying various expressions. Ideally 
one would like to be able to do this integral at the level of symbols. We hope that this 
integral expression is of use in determining the symbol. If we perform the OPE expansion, 
by expanding in powers of e~^, then it is easy to do these integrals term by term since the 
projection onto positive or negative frequencies is trivial once we do the expansion in e""^. 
Thus, this expression can be used to check that a possible guess does indeed obey the OPE 
expansion. 



5 The Hexagon Wilson loop 

A simple compact expression for the remainder function for the hexagon Wilson loop was 
found in [7], based on previous work [6, 5]. In this section we rederive the Hexagon Wil- 
son Loop expectation value at two loops without using these previous results. The main 
tools/assumptions will be that the OPE expansion holds and that the full result has a sym- 
bol. There are very strong self-consistency checks of both hypothesis so to a very great 
extent they are themselves proven. 

The hexagon has three independent cross ratios ui,U2 and U3 given by 

«^ = -2 -2 (46) 

■^i-l,i+2-^j+l,j-2 

where the Xj's are the cusps positions and x'fj = (xj — Xj)"^. We have x-j = xi etc. It is also 
convenient sometimes to use the variables 61, 62, ^3 and /i given by 



h = .l ' , ^^ = A + ^/A^ (47) 

where 2 A = — ^1 — &2 — &3 and U4 = ui etc. The two loop remainder function is a totally 

symmetric function of Mi, U2 and M3. We can also express /i and hi in terms of momentum 
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twistors^ 

bi _ (Aj_3, Aj_2, Aj, Ai+i)(Ai_4, Aj_3, Ai_2, Aj_i) ^^^^ 
fi^"^^' (Aj-3, Aj_2, Aj_i, Ai)(Ai_4, Ai_3, Aj_2, Aj+i) 

There are three OPE expansion channels which correspond to the three inequivalent 
choices of two opposite edges along which we stretch the polygon. This is conformally 
equivalent to picking two adjacent legs which we make coUinear. In each OPE limit one of 
the Mj's goes to zero while the other two remain finite. For example, in the channel where 
M2 — )■ the OPE limit can be obtained by picking a reference square as indicated in [2] (see 
appendix F.3 in [2]). In terms of the symmetries of this reference square the three cross 
ratios read [2]^ 

sinhrtanhr 1 

Ui3 = e —. r, U2 = (49) 

2 (cosh r cosh a + cos 0) cosh r 

The corresponding OPE limit is the limit where r — )■ oo with a and held fixed. 

We will now describe the logic of the derivation of the two loop Hexagon remainder 
function. The technicalities of the several steps are then explained in detail in the subsections 
which will follow. 

• We start by computing the U{1) result in a given channel (23), (24). All channels are 
equivalent so we can pick one of them, say the one where U2 0, i.e. r — )■ oo with 
(49). We find^ 

ru{i) = -Li2 (1 - Mi)+Li2 (n2)-Li2 (1 - n3)+log^ (1 - n2)-log (ui) log (^3) + — (50) 
This is the seed of the computation. It is computed in greater detail in section 5.1. 

• Next we rewrite T'fyci) as a sum over particles propagating in the flux tube. At this order 
the particles are free and their energies are quantized to integers. The SL{2) symmetry 
of the two null lines is preserved at this loop order and therefore the propagation of 
primaries and their descendants is conveniently packaged into SL{2) conformal blocks. 
The relevant conformal blocks are computed in section 5.2 and read 



J-^,p(r) = cosh-2/^(r) 2F1 



/3-^|,/3 + 2|,2/3,cosh-^(r) 



(51) 



where p is the momentum conjugated to a and /3 is the SL{2)r conformal spin of the 
free particle states. The set of primaries that can flow were described in section (4.1). 

The seed (50) is neatly decomposed in terms of exchanges of these particles as 

/dp °° 



m=— 00 



-^Here (Ai, Xj, Xk, Xi) stands for Xi A Xj A A A;. 
2 One minor change (^horc = ^thoic + tt. 

^We use ^^^^r = 5^ = a/2 where a is used in [10, 7]. This means that the two loop reminder function - 
defined as what multiphes g"^ - differs from the one in [7] - defined as what multiphes - by a factor of 4, 



19 



where 

VUp) = 4 (-ir ^ ,r'2 (53) 

are the f/(l) form factors and B{a, b) = This is a very non-trivial check of the 

OPE expansion. Formula (52) is derived in section 5.3 using a very useful differential 
operator 

D = -{dl + dl) . (54) 
It would be interesting to find a nice physical meaning for this operator. 

Having identified the primaries that are fiowing we can promote our result from one loop 
to two loops by correcting their free dimension to their quantum corrected anomalous 
dimension. More precisely the two loop result when M2 — )■ is given by 

R2I00PS = ^ log(M2) D2 + R2 , (55) 

where D2 and R2 have regular power series expansions in U2 oc e~^^. D2 is the discon- 
tinuity of i?2 around U2 = 0. We can identify D2 by dressing the conformal blocks in 
(52) by the anomalous dimensions of the corresponding primaries (39). That is^ 

/dp °° 
_g-«P^ ^ COs{m(f)) Cm{p) Im+2{P) J^m/2,p{'^) (56) 
m=—oo 

In principle we could try to directly resum this expression, e.g. by making an educated 
ansatz for the kind of functions that could appear in the result. The result is given in 
formula (86). 

A neater way of resumming (56) is by assuming that the resulting function has a symbol 
and computing this symbol. Then we find the function corresponding to this symbol 
and with the right reality properties and verify that it indeed resums (56). As argued 
in section 5.4.1 the symbol of D2 is given by 

Sym {D2) = Ml ® Xi + M3 ® X3 + (1 - U2) g) Y2 . (57) 



• The functions Xi and X3 can be found by taking the limit where ui and U3 are sent 
to infinity respectively. They will be computed in section 5.4.2. 

• The function Y2 is related to the point where we reach the radius of convergence of the 
OPE expansion. That is the point where the two null lines of the OPE expansion in- 
tersect. These are the two null lines containing edges (0:3, X4) and {xi,Xq), respectively. 
Once we know Xi and X3 we find I2 using the integrability condition which ensures 
that (57) is the symbol of a function. This is done in section 5.4.3 and concludes the 
re-sumation of (56). 

''Note that 7-m(p) = 7rn+2(p)- 
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Finally we need to find a totally symmetric function of Ui, U2 and which behaves as 

1 



R2 loops = 2 log(wi) A + Ri 



, when Ui — )■ 



(5^ 



This is discussed in section 5.5. It can be done in the most natural way using symbols. 
We assume that the two loops remainder function has a symbol. Then the simplest 
possible candidate is 

Sym[_R2 /oops] = ui® Sym[Di] + M2 Sym[£)2] + U3 ^ Sym[D3] (59) 

as explained in section 2.2. 



5.1 The one loop source, r[/(i) 

The seed for the perturbative OPE expansion in a given channel is the one loop OPE 
result (23), (24). As reviewed in section 2, it is given by 



ru(i) = log 



(60) 



1/(1) 



By computing that ratio explicitly, we find the result (50). To be more precise, we find the 
result 



r{;(i) = -Lis (1 - Ml) + Lia (M2) - Lis (1 - U3) + log^ (1 - U2) - log(ui)log(M3) + 

+ log(Mi/u3)log(l -M2) 



TT 



(61) 



The expression (50) is obtained from (61) by symmetrization under ui which amounts 

to removing the last term log (ui/u^) log (1 — M2). Since that is a symmetry of the remainder 
function, such a symmetrization does not affect the final result. In fact, this last term is killed 
by the kernel that encodes the anomalous dimensions. The reason is that it corresponds to 
essentially zero momentum modes of the particle that has no anomalous dimension at zero 
momentum (as explained in [2] these can be viewed as Goldstone modes of some generators 
in SL{2)t-). Since the technical details of the computations are also relevant for the next 
section, we will now present them in some detail. 

To compute ?"c/(i), we start by making the following choice of momentum twistors for the 
reference square 



(62) 



To identify the left and right null lines of the square we construct the dual momentum twistors 
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(b) 

Figure 3: (a) We see the hexagon. The sides 1 and 4 are the left and right sides. The vertices 
1 and 4 coincide with the reference square. In (b) we see the contours that give rise to the 
U(l) source ^'[/(i). 



With this choice, the symmetries of the reference square (62) and the SL(2) symmetry 
of the two null lines act as multiplication by the SL{4) matrices 
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(63) 



where ad — be = 1. Here, t is associated with the "Hamiltonian" generator that that moves 
points on the left and right sides of the reference square towards the bottom side, s is 
associated with the "Momentum" generator that that moves points on the bottom and top 
sides of the reference square towards the left side. The coordinate / is associated with the 
SO {2) "Angular Momentum" which parametrize rotations (boosts) in the two transverse 
directions to the reference square. Up to a global S'L(4) transformations and rescaling, the 
Hexagon momentum twistors can be brought to the form 
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(64) 



The corresponding top and bottom Pentagons are 

Top = (Aleft, A2, A3, Aright; A 

bottom J 

Bottom = (Aieft, Atop, Aright, A5, Ae) . 



(65) 
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Indeed two of the null edges of all these polygons lay on the two null lines identified above. 
For example, the dual momentum twistor /iieft = Ae A Aioft A A2 oc yUioft and so on. Note also 
that the reference square (62) has two cusps that coincide with two of the hexagon cusps. 

That is, Aright A A5 = Aright A Abottom and Aleft A A2 = Aleft A Atop- 

This choice of reference square is convenient for ^(/(i) to be finite and free of conformal 
anomalies and therefore to be a function of the conformal cross ratios Ui. On the other hand, 
r, a and in (49) are associated with the three symmetries of a different reference square. 
That reference square was considered in [2] and realizes the Ui O ^3 symmetry as a ^ —a. 
The relation between the t, s, f conformal frame and the r, a, (j) conformal frame is^ 

e* = sinhr , e"* = e*^ cothr , f = i(p (66) 

The hexagon, pentagon and square polygon Wilson loop in a U{1) theory can be found in 
[10]. Plugging these into (60), using (18) and symmetrizing in ui ^ U3, we arrive at (50). 

5.2 Conformal blocks 

We interpret ru(^i) as a sum of excitation freely propagating on the flux tube. These exci- 
tations are created at the bottom and been absorbed at the top, see figure 3(b). When we 
will go from one to two loops, both the propagation as well as the creation and absorption 
of excitation are corrected. At one loop, the flux tube is invariant under the SL{2) transfor- 
mations UsL{2,R) which act on the bottom of the square. Therefore, the two loop correction 
to the propagation of an excitation due to interaction with the flux tube does not break the 
SL(2) symmetry [3]. It commutes with it. This means that the correction to the energy of 
excitation is the same for the primaries and descendants. This is why it so convenient to 
decompose the excitations of the flux tube into irreducible representations of SL{2, R). The 
blocks of excitations forming an irreducible representation are functions of t, s and / known 
as conformal blocks. Different blocks are uniquely identified (up to overall normalization) by 
the SL{2) Casimir. To compute the hexagon conformal blocks, we first need to understand 
how t, s and / transforms when we act with the SL{2) transformation in (63) on the bot- 
tom of the hexagon. We will then express the SL{2) Casimir as a second order differential 
operator in s, t and /. 

Under the SL{2) transformation Usl{2,r), s,t and / transform as 

e^*' = (a + c){b + a e'*) , e^'+*' = a e'^' , /' = / . (67) 

Since / does not transform, the conformal blocks only depend on s and t. The transformation 
(67) does not respect the naive group multiplication. The reason is the following. The 
cross ratios are SL{A) invariant. Therefore, we can not ask how they transform under a 
specific SL{2) group element without first fixing a frame. The hexagon parameterization 
(64) is one such choice of frame. After we act with gi G SL{2) on the bottom, the frame 
changes. To rotate it back to the original frame, we act with an M{gi) G SL{A) compensating 

^Actually, in [2] we really had /here = *</'therc + i""- Here we have a slightly different definition which 
introduces some irrelevant minus signs. 
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transformation and rescale the twistors accordingly.^ If we now act with a second g2 G SL(2) 
transformation on the bottom, the resuh is not the same as acting with g2gi because g2 do 
not commute with M{gi). Instead, we should act with M{g2gi) g2 M~^{gi). In other words, 
when constructing the Casimir operator, we should then use the covariant derivatives. 

There is a shortcut however. The SL{2) Casimir can be written as 

C5L(2) = -Jl + J2 - 4 , (69) 

where for irreducible representation of spin f3 we have csl{2) = (3 {(3 — 1). Since the finite 
action of any group element is given by (67), we can read the action of Jf directly by picking a 
group element Usl(2,r) = ^2x2 x expanding the result (67) in 9 and picking the quadratic 
term. We can then sum over i = 1,2,3 to get the action of the Casimir operator. We find 
the corresponding differential equation for the conformal block' 

[e-2*52 + 2(1 - e~^'ds)dt + (1 + e-^')d^ - AI3{(3 - 1)] F(s, t) = . 

It is very convenient to switch to the r, cr, variables in terms of which the a —o symmetry 
is manifest. We find 

{dl + 2 coth(2r) dr - p'sech2(r)] e'^" J}3,p(r) = 4/3(/3 - l)e'''" J}3,p(r) . (70) 

where p is the conjugate momentum to a. The solution to this differential equation with the 
appropriate large r behavior is given by (51). It is of course p -H- — p symmetric. For future 
use let us quote the large r behavior of these conformal blocks, 

1 / + »2 1 \ , , 

■^""-'W = ^^(^ + ^^rf7 + -j^ 

We will now decompose the f/(l) contribution in terms of these conformal blocks. 



5.3 Decomposition of rf7(i) using □ r(7(i) 

We will now decompose the f/(l) part using the conformal blocks (51) derived above. As 
explained in section 4.1 we expect two infinite towers of conformal primaries. In one of the 
towers the primaries have conformal spin |m|/2 and 50(2) angular momentum m. In the 
other tower the conformal spin is |m|/2 + l for the same angular momentum m. We therefore 
expect, and will derive next, that r[/(i) = r_|_ + r_, where 



-ipa 
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^Explicitly we have 



J2 e'""^CUp)J'm/2A^) + Yl e"''"'^C'm(p)^W2+i,p(r) 



m=l 



m=0 



(72) 



M{gi) 






V 





a-\-c 
a 









a-\-c 






1 / 



^Note that contrary to the M}'^ case, where the SL{2) conformal block is simply F — log(l 
s enters the differential equation. 



-2t' 



(68) 



, here 
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and r_ is obtained from r_|_ by replacing e**^ — )■ e A simple way of proving this claim 
without decomposing the result is by using the projectors 

V± = d^ + 2 coth(2r) dr + sech2(r) dl + d^d^ T 2i) (73) 

which project out each of the two towers T>±r± = 0.^ Acting with these projectors directly 
on (50) leads to 

V+V^ run) = (74) 

which implies (72). Of course, after deriving Cm{p) and Cm{p) we can trivially confirm 
directly that (72) holds. 

To derive these structure constants it turns out to be quite useful to use the box operator 
defined in (54). It is a sort of two dimensional Laplacian in the directions (j) and a and turns 
out to simplify the r^^i) dramatically. In Fourier it amounts to multiplying the integrand in 
(72) by (p^ + m^). The reason for the remarkable simplification of r[7(i) still eludes us. We 
find 

1 + sec(0j cosh((Tj cosh(rj \cosh(crj cosh(rj / 

The Fourier transform in a yields 

□-«-|/g-(-)"--(f-'|.f-«D(^ 

For large r we find therefore 



r = 4 V / ^ _ cos(m0^ (-1)- B j^,^) 
^^^^ 271 cosh™(r) 2^'»>o (p2 + ^2) 



(2 + m)(m^ + p^) ^ 
Am(m + 1) cosh^ r 



m=0 ■ 

where we inserted back the denominator {p^ + m^) to remove the box.^ Using (71) and 
comparing this expansion with the expansion of (72) we read of the structure constants. We 
find 

/dp i„„ I v-^ cosimd)) cos((m — 2)6) \ ^ , . ^ , , , 

i E -r^ + E p.|(,„_2;' CM ^™/.,(-) (76) 
\m=l m=2 ^ ' / 

where 

It is easy to see that the decomposition (76) is of the expected two towers form (72). It is 
also equivalent to the expansion (52) presented above. To show this one simply needs to use 
the identity C^ij)) J^m/2,p{j) = C2-m{p) ^2-m/2,p{j) whicg holds for integer m. 



^Also 'D±rzf: = ^Aidcf, rip. The building blocks r± should be identified with r+ and r_ discussed in section 
3. See in particular the discussion in section 3.5 which leads to this identification up to simple terms in the 
Kernels of the differential operators 'D±. 

^At the end of the day we directly check the decomposition and thus verify that we did not loose any 
zero modes. 
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The structure constants Cm{p) are meromorphic functions which vanish very fast at in- 
finity and which have infinitely many poles along the imaginary axis. We can compute the 
integral (76) by residues and hence obtain, for each m, an expansion in powers of e"'*^' as 
expected. 

5.4 Two Loops Resummation 

Now we will promote our one loop seed to two loops. When going from one two two loops 
three types of corrections. The particles energies acquire anomalous dimensions, the form 
factors are corrected and finally we have two particle exchange. The energies appear in the 
exponent multiplying r oc log(M2). Therefore, when we expand the result in powers of the 
coupling, the anomalous dimensions of the excitations, which appear in the exponent, will 
lead to a contribution to the remainder function which is linear in r. The other two effects 
do not contribute to the r linear part. That is, at large U2 we have 

R2I00PS = ^ log(M2) D2 + R2 , (78) 

where D2 and R2 have regular Taylor expansions in powers of U2- Using the OPE expansion 
we can compute D2 by simply dressing the expansion (76) by the anomalous dimensions 
(39) (40) of each of the conformal blocks, 

^ [dp i ( ^ 7m+2(p) cos(m0) ^ lra-2{p) cos((m-2)0) ^ 

^' = J 2^' p- + m^ + p2 + (m - 2)2 ) ^^^^^ ^^^''^'^^^ ' 

(79) 

The challenge is to perform this sum. 

There is a simple brute force way of doing it if we assume that the result is an expression 
of transcendentality three given by a bunch of logarithms and polylogarithms. There are 
only a finite number of reasonable arguments which we can expect to encounter since we do 
not want to have discontinuities or other singularities at non-physical loci. Hence any such 
ansatz will have a very large but finite number of terms. We can then expand the ansatz at 
large r and match it with the expansion of (79) to fix the coefficient of each of the logarithms 
and polylogarithms. This would be the analogue of the unitarity method approach where 
one makes an ansatz for the result using a basis of box integrals whose coefficients are fixed 
by matching with the expected discontinuities of the result. 

We will follow another approach which uses the technology of symbols to break the 
computation of the sum (79) into much simpler blocks. 

5.4.1 The symbol of D2 

We assume that D2 has a symbol. At the end of the day, when we will compute D2 and 
compare directly with (79) we will check that this assumption is correct but for now let us 
motivate it and understand its consequences. 
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Under the assumption that D2 has a symbol, that symbol should be of the form (57) 

Sym[D2] =ui® Sym[Xi] + M3 ® Sym[X3] + (1 - U2) ® SymiFa] (80) 
as we now argue. 

First note that the function D2 is the discontinuity of the two loops remainder function 
around U2 = 0, where edges 1 and 2 become coUinear^*^. What singularities/discontinuaties 
can D2 have? First, it cannot have another discontinuity at U2 = 0. Such a second disconti- 
nuity would mean that has a term linear in r in contradiction to the OPE expansion. 
D2 can have a discontinuity when two other edges become collinear. These are the points 
where edge 2 and edge 3 become collinear^ ^ or when edge 3 and edge 4 become collinear^^. 
At these points M3 — )■ and Mi — )■ correspondingly. 

The other singularity that D2 may have is at the radius of convergence of the OPE 
expansion. That is the point where the two null lines of the OPE expansion intersect. These 
are the two null lines containing edges (x3,X4) and {xi,xq), respectively. At that point 
U2 = I. This self-crossing does not happen in the Euclidean sheet. These three possible 
singularities are precisely the ones captured by (80). 

The only other singularities/discontinuaties one may expect are when two cusps Xj and 
Xi+3 become null separated. At these points, one of the Uj's diverge. These are however 
already included in (80). 

We conclude that to derive (the symbol of) D2 we can focus on finding (the symbol of) 
Xi, X3 and Y2. This is what we do in the next two subsections. 

5.4.2 The symbol of 

In this section we focus on computing the contributions Xi and X3 to the sum (79). They 
are related by symmetry, Xi = X3|^^___^^, and therefore it is enough to focus on X3 which is 
what we will now do. 

The contribution to X3 comes from the large a behavior of D2. More precisely, we should 
compute the a linear term of D2 when o" — > 00. 

Those come from double poles in p. The residue of a double pole is the derivative of what 
multiplies the double pole. This derivative will, in particular, act on the exponential e"*^'^ 
bringing down an —ia factor. We are not interested in the the action of the derivative on 
the rest since it will not lead to a contribution linear in a. That is, we can simply replace 
the double poles by —ia times the factor which multiplies the double pole evaluated at the 
position of the double pole. Then we sum over all possible double poles. In this way we 
compute the term linear in logM3. The details of the summation are given in appendix B.l. 

^''or, conformally equivalently, 4 and 5 
^^or, conformally equivalently, 5 and 6 
^^or, conformally equivalently, 1 and 6 
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We find 



Xs = -2U2I — ^ ) -2Li2(-f ) -Li2(l-ni)- Li2 (1-1x2) -Li2(l-M3) 



fibi J \ bi 

-log^(Ml) -log(l -U2)log ( — ) -log(l -M3)log ^2) U3 



Ui J \ Ui 

Xi is obtained from X3 by ui ^ U3, bi o 63. At the level of symbols we have 

-Sym Ai = 2— ® ® h ^3 ® 

2 03/^ + 1 M2 U3 ui U3 

-SymXg = 2—^ 7^®/^ ® ® h Mi ® 

2 Oi /i + 1 Ml M3 Ui U2 Ui 

5.4.3 The symbol of Yi 

The goal of this section is to constrain the remaining building block, Y2. Since Y2 is the 
discontinuity of the result at the boundary of the convergence region it is not straightforward 
to compute it using the expansion (79). Instead we will finding using some symbol technology. 

To start we need an ansatz for Y2. Looking at the symbols of Xi and X3 we see that 
in their last slot there is always either fi ot 1 — 1/uj. Recall that derivatives of the symbol 
act on the last slot. Zeros and poles in the last slot translate into pole singularities for the 
derivative of the corresponding function and therefore have a physical meaning. The simplest 
ansatz for Y2 which does not introduce new singularities is therefore 



1 ^ 

-Sym[r2] = /o ® /i + 5Z /. ® (1 + l/u,: 



i=l 

We will see that this ansatz is indeed general enough. To determine /o,i,2,3 we impose that 
Y2 is a symbol of a function (10). That condition almost fix Y2 completely. The details are 
given in appendix B.2. We find that if and only if 

1^ .1 1 — Ui 1 — M3 , X 1 — U2 , s 

-Sym[y2] = -Ml ® M3 ® - 4^ ® + /(M2) ® (82) 

2 Ml M3 M2 

then Sym[y2] is the symbol of a function Y2. The unfixed function f2{x) should be a rational 
function which we will now fix. 

To find f2{x) we use the information that D2 is given by a sum of the two towers of 
conformal blocks. In other words, 

V_V+D2 = (83) 

where V± project out each of the two towers, see (73). The action of V_V^D2 yields a term 
proportional to log(l — M2) plus a term without any log. Both should be zero separately and 
hence we obtain two differential equations for f2- One of them is 

Q = 2xg" + g'-^^^—^^ g = \ogf2 (84) 
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The other is a third order differential equation^^. The solution of (84) is 

/2(^)=(r^C7,e-^^/v^ (85) 

The constant C2 is irrelevant for the symbol and we must set Ci = since /2 should not 
be a transcendental function. Indeed, plugging (85) in the above mentioned third order 
differential equation we see that we get zero iff Ci = 0. Therefore, we obtain 

1^ r -, 1 — Ui , , 1 — Uo 1 — 1 — Ms 

-Sym I2 = -ui ® 2 (1 - M2) ® + M2 ® M3 ® 4/i (g) /i . 

2 ui U2 U2 U3 

5.4.4 Final expression for the discontinuity D2 

Once we have the symbol of the discontinuity D2 it is easy to construct a function whose 
symbol is the same as the one of D2. The function D2 is then determined (up to a constant) 
by having only physical branch cuts [7]. We find^'* 



Ml 
U3 



2 log 

log (1 - U2 



Li2 



+ Li3 

1 

fibi 



/i 



Li2 



- 2 



Lis ( -— ) + Lis (-/W&2) 



+ 



2Li2 1 



Ui 



- log (M1M3) 

1 



Lio 



+ U2 (-Ai&2) 



Li2 {U2) + - log (1 - M2) + 2 log 



1 2 f ^ 



Us 



+ 



+ 



2 log (J- - 1^ log (Ml) log (M3) I + { 1 ^ 3} 



(86) 



It is indeed straightforward to check with Mathematica that this function indeed resums 
(79).^'' Equation (86) coincides indeed with the discontinuity at M2 = of the result in [6, 7]. 



5.5 Symmetrization and full result 

So far, we have computed D2, e.i. the term linear in r of the two loop remainder function 
at large r. In other words, we computed its discontinuity in M2. The two loop remainder 
function we are interested in R2, is a symmetric function of the conformal cross ratios mi^2,3 
(that is cyclic and parity invariant) hows OPE expansion at small M2 is given by D2 (86). That 
information is not sufficient to uniquely determine i?2- The situation is analogous to trying 

i3lt is given by = g'" + i0t)9" + Vil-l)^ 

^^The conformal cross ratios do not fix the polygon uniquely. As a result, the remainder function and, 
in particular, D2 are multi-valued functions. Whenever we write such multi-valued function explicitly, we 
mean in the Euclidean sheet, namely, when all distances, that are not automatically null, are spacelike. 

""^^This can done either numerically with very high precision or by expanding this functions to extremely 
high orders in l/cosh(r) and comparing with the explicit expansion (79). 
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to determine a four point correlation function at some loops level from its discontinuities. 
On the other hand, if we assume that R2 has a symbol, then there is a natural candidate for 
its symbol 



Sym[_R2 ioops] = ui® Sym[Di] + U2 ® Sym[£)2] + ^ Sym[D3 



^7) 



The structure of the symbol (88) is manifestly symmetric in the Wj's and was motivated in 
section 2. 

More explicitly, using (80) we see that the structure of the result is quite neat. 



Sym[i?: 



2 loops} 



Ui Uj (g) Xij + ^ (g) (1 - Ui) (g) Yi 



In the previous subsections we computed X21 and X23 (which we denoted as Xi and X3 
respectively) and Y2. All other Xij and Yi are related to those by trivial relabelings. The 
several OPE limits commute and this implies X^j = Xji which is indeed a property that we 
can directly check from the solutions that we get. 

A simple but important self-consistency check is that (88) is a symbol of a function. A 
priory that is not the case for generic Di. To check that (88) is a symbol of a function note 
first that since the functions Di (86) do have a symbol (80), it is enough to apply the check 
(10) to the first two slots of (88). Applied to the first two slots, the check (10) yields zero 
provided Xij = Xji. That is indeed the case and follows from the the fact that the different 
OPE limits commute as we just explained. 

Again, up to a constant, we can find the function whose symbol is given by (88) by 
requiring the presence of physical branch cuts only [7]. The resulting two loops remainder 
function is then given by [6, 7] 
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19) 



6 Conclusions, discussion and speculations 

In this paper we have computed the discontinuities of the two loop remainder function for 
polygonal null Wilson loops for a general number of sides. Here we considered the "OPE 
discontinuities". These are defined as follows. We pick two sides of the null polygon, which 
divides the polygon in two halves. We act on one half of the polygon with the operator 
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e'^^ where H is a "Hamiltonian" in SL(4) which preserves the two null sides, then this 
Hamiltonian squeezes this half of the polygon. We can perform the full OPE expansion of 
the one loop answer, keeping all powers of e""^. This is what we call the one loop result or 
rf/(i), since it can be computed in a free U{1) theory. The OPE expansion is defined for any 
loop order. So for any loop order we can take a contour in the space of cross ratios where we 
deform the contour by setting a large value of r, then we take r — )■ r + 27ri. In a perturbative 
expansion, this picks out all the logarithmic terms that arise due to anomalous dimensions. 
We call this "the OPE discontinuity". In the case of the two loop answer, there is only a 
linear term in r. The coefficient of this linear term in r is the discontinuity D of the two 
loop answer. It is a transcendentality three function. It can be computed by considering 
the one loop OPE and then multiplying each term by the anomalous dimension. Since the 
anomalous dimension depends continuously on a "momentum" quantum number, we find a 
final expression for D which is has the form of a convolution of a simple kernel acting on the 
one loop expression. The final answer, given in (43) (45), is relatively simple and compact. 
Note that in this form, it is quite manifest that the final answer is a transcendentality three 
function. The one loop answer is a transcendentality two function and the integral (43) raises 
the transcendentality by one. The integral that projects onto positive or negative frequencies 
(44) does not raise the transcendentality because it is a closed contour integral. Note that 
the transcendentality of the answer is not manifest in any of the previously known integral 
representations for the Wilson loop or the amplitude. 

It would be nice to be able to evaluate the symbol of (43) or (45) without having to 
evaluate the integrals. A direct evaluation of the integrals seems possible and doable (as 
long as one has enough patience with mathematica...). On a term by term basis, these 
integrals seem to give square roots of momentum twistor cross ratios. We do not know if 
they all cancel in the final answer, since we have not computed all the terms. 

We have argued that the knowledge of all OPE discontinuities should be enough to 
reconstruct the symbol of the two loop answer, for any number of sides. 

There are some terms in the symbol of the final two loop answer which we can predict in 
a simple way. Namely, we can have terms of the form Uij ® Uki ® Xij^ki + Wij ® (1 — Uij) ® Yij 
where i, k,j, I label sides of the polygon which are cyclically ordered. Here Uij are the cross 
ratios made with the four points at the ends of the two null sides. These terms arise from 
"crossed" propagators. These are terms which appear in the remainder function due to an 
oversubtraction when we say that W = e'"™''p"^i-'°°p+^. These are terms that are present in 
the exponential of the one loop answer but are not present in the planar two loop answer. 
Thus they contribute to the remainder function. For some such terms, we can even say what 
Xjj fc/ and Yij are. If i,j,k,l are well separated from each other, these crossed propagators 
are the only contribution to the symbol. Now, in the case of the hexagon, there are other 
Feynamn diagram contributions. However, we have experimentally seen that the symbol of 
the two loop answer still has only these two type of terms in the first slot. Some of us want 
to conjecture that this is true in general. Other of us want to be more cautious since the 
evaluation of individual terms in (45) give terms which are not of this form (which could 
eventually cancel). Clearly, if one makes a correct hypothesis regarding the possible terms 
appearing in the symbol, then the computation of the answer reduces to a finite problem, 
determining the precise combination of a finite number of possible symbols. In that case. 
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the explicit OPE discontinuities we have discussed could give a large number of constraints 
which would help in determining the answer. 

We have shown how to derive the answer for the hexagon using our method. In fact, in the 
hexagon method we did not follow word by word the general route outlined in the first couple 
of sections, which we could have also done. Instead, we performed some explicit Fourier 
transformations, which lead to particularly simple formulas. In particular, we noticed that 
the action of a particular differential operator (□), simplified the expressions considerably. 
It would be nice to see if this simplification extends to other cases or whether it is just an 
accident for the hexagon. It is reminiscent of the recursion relations for integrals studied in 
[12]. In the hexagon case, one could fix the discontinuity by looking at its own discontinuities. 
Of course, one would like to be able to do this for the general case, so as to compute the 
symbol step by step. 

We also considered the heptagon Wilson loop (in appendix A. 2). We have derived its 
OPE discontinuity in a particular mixed kinematical regime. That result can be used in the 
future to restrict the form of the heptagon two loops reminder function. 

As we argued in [3], we expect our two loops results to be valid for null polygon Wilson 
loops in any conformal gauge theory with a weak coupling limit. 

Our expression for the discontinuity amounts to an insertion of the anomalous dimension. 
In fact, we could imagine we inserting e^'^*^ — 1, where H is the full Hamiltonian of the 
interacting theory. This gives us the discontinuity of the answer. It also has the interpretation 
of adding (a combination) of conserved charges along the worldsheet (or planar diagram). 
From the point of view of the amplitude, these are non-trivial integrability charges. It should 
be interesting to understand how to add more general integrability charges in order to further 
constrain the problem. 

In the particular case of three loops, we can find some of the discontinuities in the same 
way. For example, in the OPE limit, there will be terms which are given by applying 
the anomalous dimension kernel twice to the one loop answer. This constrains some of the 
terms in the symbol. In particular, we can determine terms like ® ® S, by 

considering the OPE that leads to a collinear limit. In addition, we can do the following. 
Imagining that one had solved the two loop problem completely. Then one could find the 
non-logarithmic parts of the two loop OPE. This involves the exchange of two particle states. 
Then one can multiply these by the one loop anomalous dimensions and, together with the 
two loop anomalous dimension acting on the one loop term, we can get the linear in r terms 
in the OPE. 
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A Simple examples 

In these section we will demonstrate by explicit examples how our techniques can be used 
for obtaining new results at two loops. We will consider two examples: The first is a double 
scaling limit of the hexagon where a single primary per 5*0(2) angular momentum propagates 
in the flux tube. The other is the full OPE discontinuity of the heptagon two loops remainder 
function in a mixed kinematical regime. Both examples can easily be extended to any number 
of edges. 



A.l Double scaling limit 

In this section we will consider an interesting kinematical limit where the two loop com- 
putation simplifies dramatically. Remarkably the final form of the result in this limit 
captures most of the features of the full result. The regime we want to consider is that when 
T — )■ +00 diverge with ( = t ~ held fixed and large. In this limit ui and U3 are finite 
while M2 — )■ so that the U{1) result for the hexagon becomes 

^2 



— - log(Mi) log(M3) - Li2(l - Ml) - Li2(l - U3) 




(90) 



In this double scaling limit the second tower (i.e. sum) in (72) and (76) can be dropped 
and in the first sum we keep only the terms with positive S0{2) charge. We can replace the 
conformal block by its leading large r asymptotics. For example, 



J 271 ^-^ + 

m=l 



(91) 



Physically, in that limit only primaries survive and, furthermore, only the primaries of one of 
the towers. That is, for each S0{2) angular momentum we have a single primary propagating. 

Next we want to promote this result to two loops by multiplying each summand by 
7m+2(p)- In this way we compute the double scaling limit of D2. Alternatively, in the a- 
space we can apply the anomalous dimension kernel as a convolution as described in section 
4.2. In the scaling limit only one tower with e*™""^ with m > is propagating. Therefore 
we do not need to project the result into positive and negative frequencies and the result of 
applying the anomalous dimension kernel is very simple 
00 

D2 {a, C) = ^g^ / [2 ru(i) (^, C) - run) {cr~t,C + t)- fun) (^ + t, C + t)] (92) 
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or rather of a particular symmetrization of the result in this limit 
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It is straightforward to plug (90) inside this integral and perform the integral using Mathematica 
for example.^' We find 



Li3(-i)-Li3 (-^+Li3 (-i 



2Li2 ( -62 ) log (M1M3 



(94) 



+ 2 



Lig 



1 



Li2 - 



1 



log( — ) +2l0g[(l+^)(l + ^)]l0gUil0gM3 

U3 / bi 63 



where bi = 63 = b^/fi and 62 = ^2/^- We also used the equation (I + 62) = + + 
1/63), which is true in this limit. The double scaling limit is a neat example when the one 
loop result can be promoted to a two loop prediction by means of a single and very simple 
integral which increases the transcendentalities of the functions in the U{1) result by one 
unit. 

It is curious to note that the double scaling result roughly captures half of the full result. 
It would be interesting to understand if it can be promoted in a simple way. This might be a 
shortcut towards the full result which might avoid the re-summations described in the main 
text. More explicitly, note that if we start from (94), add the same but with — )■ we 
are missing the third line in (86). 

The scaling limit can be easily generalized to other polygons with more edges hence 
providing very non-trivial constraints on the form of the discontinuity of two loop MHV 
amplitudes for more than six gluons. A simple way to state the generalization is to say 
that we choose a special polygon where all the momentum cotwistors (say Aj) of the top 
part of the polygon are the same and, in addition, they are orthogonal to A^^, \r which 
are the momentum twistors of the two special lines that we are using to extract the OPE 
discontinuity. Note that it makes sense to take this special polygon only after extracting the 
original OPE discontinuity (the linear term in r in the OPE). For such polygons, only r+ is 
nonvanishing, and r+ only has positive frequencies. Thus we can simply do the integral in 
the first line of (43) to get the discontinuity. 



A. 2 The Hybrid 

In this second example we consider the OPE discontinuity of a family of polygons we call 
Hybrid's. An Hybrid is a polygon where the top part is in M}'^ kinematics while the bottom 
is in general M}'^ kinematics. The physical reason for a big simplification in that regime is 

that the top part can only absorb type excitation and therefore we will not have two 

infinite towers of conformal blocks. Still, the result is not trivial and, in particular, involve 
polylogarithms and not just logarithms as in strict M^'^ kinematics [13, 4, 3]. 

^^Thc kind of integrals we need to do are 

/ /"^ ,^ log(aw + 6) log(cu; + d) and / ,. Li2 f 1 + (93) 

Ji 'w(w — l) Ji w{w — l) \ b + wj 

which are quite simple to handle. 



34 




Figure 4: A polygon in a mixed kinematics called Hybrid. In the OPE limit u ^ oo the 
two loops remainder function can be divided into two pieces as R2 = t D + D. The function 
D{s,t) is computed in this section. 

The simplest possible Hybrid is an heptagon where the top part (in M^'^) has three edges 
while the bottom part (in M}'^) has two edges and is the same as the bottom part of the 
hexagon (see figure 4). 

This Hybrid has only two independent cross ratios. One way to see that is as follows. 
Given an heptagon hybrid, we can use conformal transformations to fix all points on the top 
as well as the bottom point which is shared by the reference square. To fix the remaining 
bottom point which is on one of the two null lines (the cusp A7 A Ai) we need one parameter. 
We are left with the point outside the plane, in the bottom (that is, the cusp Ag A A7). It is 
null separated from the two neighboring cusps and therefore depend on two more parameters. 
However we can use S0{2) orthogonal transformations - which leave all other points invariant 
- to kill one of these parameters. Hence we have 1 + 1 = 2 independent cross ratios. These 
two cross ratios read 

y _ ^2s _ (Ai, Ae, A2, A3)(A2, A7, A5, Ae) ^ _ ^2* _ (Ai, Ae, A4, A5) (Ai, A2, A5, A7) ^^^^ 
(A5, Ag, A2, A3)(A6, A7, Ai, A2) (A7, Ai, A5, A6)(Ai, A2, A4, A5) 

The two loop remainder function is given by 

R2{s,t) =tD{s,t) + D{s,t) (96) 

where both D and D have regular expansions in l/?7 = e~^*. In this section we predict the 
function D{s,t). The large t limit corresponds to the channel in figure 4 . In this channel 
the rc/(i) is 

run) = 2Li2(-l/x) - 2Li2(-l/2/) (97) 

where 

rj — 1 

We have^^ 

1 f Q-'i-P^ — g-«ps 
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Integration slightly below the real axis. 
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Note that the role of the second term in (97) is minimal in the OPE expansion: When rj ^ oo 
we have x y and therefore the second term is only canceling the rj independent constant 
coming from the OPE expansion of the first term. Inside (99) we recognize the F+_ form 
factor [3]. 

We will now argue that the term linear in t of the two loop result has the remarkably 
simple form 

D = -F{x) + F{y) (100) 

where 

F{x) = Li3(-x) + 2Li3(^)-Li,(-x)log(^^^) (101) 

+ -logf ^^jlog'(x + l) + ylog(a; + l) 

As before the result (100) is given by a difference of a function of x and a function of y. As 
in the U{1) case, the function of y is simply regulating the result as ?7 — )■ oo and therefore 
X ^ y. 

To learn what type of excitation are propagating, we act on with the SL{2) Casimir 
operator. The corresponding differential equation reads 

VJ^^{r],y) = [r/(r7 - l)d'^ + (2r/ - l)d, + yd.dy] J^f,{v,y) = P{P - l)Mv,y) (102) 

If we act with the differential operator V on the first term in (97) we get zero. This is 
consistent with our physical picture: It means that (100) is describing the exchange of a 
single primary with (3 = 1 plus its descendants. This primary is nothing by ! As argued 
above, the second term in (97) is simply regulating the OPE expansion and therefore we will 
not care about it; in particular it does not carry any t] dependence. 

Having identified the primary which is flowing we can now promote the result to two 
loops. We simply need to multiply the Fourier decomposition (99) by 72 (p) (39) and perform 
the integration. Equivalently, we can use the convolution kernel for 72 (p). Doing so, we 
find (100). More precisely, the OPE expansion tells us to dress the first exponent in (99) by 
72 (p); this leads to F{x). We can think of —F{y) as regulating the OPE limit rj —¥ 00 which 
leads to X ^ y. 

The computation of the full remainder function for the Hybrid does not seem, however, 
to be considerably simpler than the computation of the general heptagon. The reason is 
simple; the OPE expansion in other channels will take us out of the hybrid kinematics. In 
those channels we will have all primaries propagating and not just 

B Technical details for section 5 
B.l More details on 

In this section we compute the discontinuity of D2 around ^3 = 0. That is the function X3 
defined as the a linear term of D2 when a — > +00 (80). This term comes from double poles 
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in p of the Fourier transform of D2 (79). The residue of a double pole is the derivative of 
what multiplies the double pole. This derivative will in particular act on the exponential 
g-ipcr ]3ri]2ging down an —ia factor. We are not interested in the the action of the derivative 
on the rest since it will not lead to a contribution linear in a. That is, we can simply replace 
the double poles by —ia times the factor which multiplies the double pole evaluated at the 
position of the double pole. Then we sum over all possible double poles. In this way we 
compute the term linear in a. 

For each m there are infinitely many poles. More precisely, for each m, we have the 
simple poles in the structure constants Cm{p) - which are already present in the ^[/(i) result 
- multiplying the simple poles in the anomalous dimensions 7m(p) 

Ai 

p — i[m + 2k) 

These give rise to infinitely many double poles at integer imaginary values. These are not 
all double poles but almost. The remaining ones are very simple and finite (for each m) and 
will be identified and discussed below. 

The key observation is that the poles of the anomalous dimension always have the same 
residue since they are made of polygamma function which are log derivatives. As seen in 
(103) the residue is 4z. Hence, from the discussion above we conclude that the effect of the 
'~fm{p) functions is trivial and simply gives 4z times —ia times the result without the gamma 
functions which is nothing but the U{1) result! This is 



;-2)X3 = -4rc/(i) + 



contribution from finitely many 
poles which we will now discuss 



(104) 



Let us now turn to the poles which we did not consider. They come from two places: 
1. The pole at p = im for m > 1 in the first tower in (79) 

T^. . . [dp -ipa^ 7m+2(p)cos(m0) 

First tower = —e^ } — Crn{p) J^m/2,p{'^) ■ (105) 

m=l 



The simple pole of Cmip) at p = im is not dressed by a pole of ■jm+2{p), so we overshoot 
when writing ^[/(i) in (104) since ?";7(i) contains all poles. On the other hand l/(p^+m^) 
does have a simple pole and therefore there is indeed a double pole at j9 = im in (105). 
So we need to sum the contribution from the double pole from l/(p^ + m^) and subtract 
the U{1) contribution from the double pole at p = im. The resulting contribution is 

= 8L. + 8L. f-ii V 4 log ( <'-^""--^> ) lor ^ - ^> 



bin J \ bi J \ Ui J \ u\ 

and should be subtracted in X3 (104). 

The double pole at p = for m = 2 in the second tower of D2 (79), 

Second towerU^2 = / ^ ^ '^^(p) J^iA^) ■ (106) 
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Here, we have a double pole at p = in r(/(i) but no pole at all for the dressed 
result D2 since the anomalous dimension vanishes quadratically in this case. Therefore 
we do not pick that contribution from the U{1) result. The contribution is given by 
£2 = -2aC2(0)J-i,o(r), 

£2 = A\og{^^ log(l-n2) . 

It is independent of (p and should be subtracted from r[/(i) in X3 (104). In other words 
we should remove the independent a linear contribution in rui^iy 



All together, we get that {—2)X^ = — 4rc/(i) + Si + S2, i-e. 



X, 



Ui 



-21og2 (Ml) -21og(l-M2)log 



U2 
Ui 



21og(l -M3)log 



Ml / 3 



B.2 More details on Y2 

As argued in section 5.4.3, Y2 takes the form 

3 

Sym[r2] = 2/0 ® ;u + 2 J] /, ® (1 + 1/m,) . (107) 

i=l 

We will now proceed and find the functions /o, /i, /2, /s by imposing is that Y2 is the symbol 
of a function. To check that Y2 is a symbol of a function we need to pick a set of independent 
variables. A very useful set of such variables is 

XjEiiiMM)- (108) 

The Ml ^ M3 symmetry of Y2 leads to /i(/i,&i,&3) = /3(/i, 63, 61). The checks that Y2 = 
a®h®c + a'®h'®d... is the symbol of a function function read (??) 

= {dj log a dk log h — dk log a dj log b) log c + {dj log a dk log b' — dk log a dj log b') log c' + . . . 
= log a {dj log 6 dk log c — dk log 6 9^ log c) + log a' {dj log 6' 9^ log c' — dk log b' dj log c') + . . . 

Each of these equations gives us several constraints since the coefficient of the different 
logarithms in the right hand side of these equations must be zero seperately. More explicitly 
we find 

= (. . . ) log 61 + (. . . ) log(l + 61/i) + . . . (109) 

for each equation. In total we find a set of partial differential equations for the functions fa- 
For example, the coefficient of log(l + bifi) of (109) with j, k = 1,2 yields 

d 

0={fi^- 1) fi {bi, 63, /i) - /i + &3/i + Ai' + 1) ^/i {bi, 63, /i) (110) 
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which imphes that 



63, //) = ( 61 + 63 + + - ) 63) 



(111) 



Wc pick this solution, plug it into the remaining differential equations, Simplify them, pick 
the simplest differential equation, DeSolve it, plug the solution again, etc. At the end of the 
day, the eight differential equations that we found almost fix Y2 completely! We find 



Sym[y2] = -2 



Ml 



Ml 



1 — ^3 \ 1 

+ M3 ® — h 4// (8) // - J{U2) C 



U2 



U3 



U2 



(112) 



where I2 is a symbol of a function for any function /2(m2)- This function is fixed in the main 
text. 



C Limiting values of r(X, Y) 



The pairing rjj^i) (X, Y) does not diverge if a single vertex of X and a single vertex of Y 
become null-separated, i.e. Xa • Fc = while all other space-time and twistor products are 
finite and non-zero. Indeed, the discontinuity A^c is zero as • Fc = 0. We can easily take 
the limit Xa -Yc ^ Q ol the answer. It takes a little bit of work if we use the expression in 
terms of space-time cross-ratios. Indeed the argument of dilogarithms where Xa ■ Y^ is in the 
numerator of the cross-ratio go to 1 in the limit, and we can drop them, but the argument 
of dilogarithms where Xa ■ Yc is in the denominator of the cross-ratio goes to infinity, and we 
need to manipulate them further. Up to constants, we can rewrite 



-Li2 



{^i ■ Ffe)(Xj_i ■ Ifc+i) 



as 



Li2 



1 - 



1 
2 



+ 7: log' 



i^i ■ Yk){Xi_i ■ Yfc+i) 
i^i ■ Yk+i){Xi_i ■ Yk) 



(113) 



(114) 



in order to bring Xa • Yc in the numerator of the cross-ratio in the dilogarithm, which can 
then be dropped. The squared logarithm terms collect to 



1, 2 (-^»-n-l)(^a-l-ig 

2 {Xa■Yc){Xa-l■Yc^^) 



2 {Xa+l ■Yc+i){Xa -Yc) 



and must be combined with the terms in r2 which contain log Xa • Y^ 



-logXa-y^log 



i^a ■ Fc) {Xa+1 ■ Yc 



C+ly 



{^a ■ ^c+l)(-^a+l ■ ^c. 



-iogX„_i -yc-iio^ 



{^a-l ■ Fc_i)(Xa • Yc) 
i^a-l ■ ^c)(^a • ^c-l) 



and 



- log Xa ■ Yc-i log — — — r - log Xa-1 ■ Yc log 



{Xa ■ Yc){Xa-\-l ■ Yc-i) 



i^a-l ■ Yc+i){Xa ■ Yc 



(115) 



(116) 



(117) 



The \og^ Xa ■ Yc terms cancel. So do the terms proportional to logX^ • Yc. What is left is a 
finite combination of logarithms. 
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On the other hand, if we use the expressions, the work is done for us: indeed none of 
the twistor inner products goes to zero, all that happens in the limit is that the argument 
of some dilogarithms goes to 1. 

We will encounter a more intricate limit, where a vertex of one polygon is brought to 
the null line defined by two vertices of the other polygon, see figure 1(b). Say that Xa is 
brought to the line which passes through and + 1. This is a more subtle limit, as 
many quantities go to zero at the same time: Xa ■ Y^ and Xa ■ ^c+i, but also (Aa,/ic+i) and 
(Aa+i, /ic+i), and the dual products. 

Let's consider the limit of r^{X, Y). It is useful to express the limit as the combination 
of three scalar constraints. Two are the obvious (Aa,/ic+i) — ^ and {\a+i,fic+i) 0. They 
are equivalent to the Xa ■ fuc+i = constraint. A third constraint follows from Xa A ^c+i = 
(notice the lack of hat!), and takes the form Xa ■ {ftc A fic+2) — ^ 0, or (Aa, /ic)(Aa+i, /ic+2) — ^ 
(Aa+i, /ic)(Aa, /ic+2)- Gcnerically, none of these four products go to zero, and hence their 
cross-ratio will go to 1. 

(Aa, Ac)(Aa+l, Ac+2) 



(Aa+l, Ac )(A 

It is interesting to look at the simplification of the following crucial ratio 

(Aa, /ic+l)(Aa+l, /ic+2) — fJ'C+l){^a, fl'C+2) 



t 



Xa ■ Yr 



c+1 



Xa-Yc i^a, P'c){^a+1, P'c+l) — (Aa, /^c+l) (Aa+l , Ac) 

Due to the third constraint, we have 

(Aa+l, Ac+2) 



t 



(Aq+I, Ac 



t 



'c+2) 



(Aa, A^ 
(Aa, Ac) 



118) 



(119) 



(120) 



There is a simple geometric meaning to the ratio t which enters the limit. The point Xa in 
the limit can be written as Fc+i ~ tY^. 

If we look at the explicit expression for r+, we have some dilogarithms whose argument 
goes to zero 



-Li-) 



(Aa+l, Ac+l)(Aa+2, Ac 



Li-) 



(Aa-l, Ac+2)(Aa, Ac+l) 



.(Aa+l, Ac )(A a+2 , /^c+1 

Some other depend on the ratio (Aa+i, Ac+i)/(Aa, Ac+i 



. (Aa-l , Ac+l) (Aa, Ac+2) . 



fl211 



-Li-2 



(Aa, Ac+l) (Aa+l, Ac> 
(Aa, Ac) (Aa+l, Ac+l) 



Li-2 



(Aa, Ac+2) (Aa+l, A- 



C+ly 



(Aa, Ac+l) (Aa+l, Ac+2) 



fl22) 



but due to the third constraint, their product of their arguments go to 1, and we can replace 
them with 



(Aa, Ac+l)(Aa+l, Ac 
(Aa, Ac) (Aa+l, Ac+l). 

Finally, we have some dilogarithms whose argument diverge. 



-Li-) 



(Aa-l, Ac+l)(Aa, Ac) 
.(Aa-l, Ac )(A a. 



Li') 



(Aa+l, Ac+2)(Aa+2, Ac+l) 
(Aa+l, Ac+l) (Aa+2, Ac+2). 



(123) 



(124) 
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and need to be manipulated by the replacement Li2{l/x) — > —Li2{x) — | log^ x, giving 



(Aa-l, /ic+l)(Aa! Ac) 
(Aa-1, /'■c)(Aa, Ac+l) 



(Ag+l, /ic+2)(Aa+2, Ac+l) 
(Aa+1, /ic+l)(Aa+2, Ac+2) 



(125) 



Now we need to combine all those squared logarithms with r^. Thanks again to the 
third constraint, all terms containing either log{Xa, fic+i) or log{Xa+i, fic+i) cancel out, and 
the result has a finite limit. 



Hence the correct limit for is obtained by adding together the main pieces 



i.k 



' = yilog(Ai,/ifc)log— ^ 



(A.,. /"/A:+i)(A,:^i. /ffc) 
(Ai, /ife)(Ai+i, /ifc+i)_ 

(Ai, fik+i){Xi+i, fik) 

P'k){\+1, P'k+l) 



(126) 



(127) 



where each of the sum omits the six terms where i — a — lori — aori — a + 1 and k — c 
or k — c + 1, and the remaining logs 



rj = log(Aa_i, i^c+i) log 7^ — r ^7^ ^ + log{Xa+2, l^c+i) log 



(Afl, Ac+2)(Aa-l, Ac) 



(Aa+2, Ac+2)(Aa+l, Ac 



and 



-^log^ 



(Aa-l,Ac) 1, 2 (AaH-l,Ac) , 1, 2(Aa+l,Ac+2) 1, 2 (Aa+2,Ac) 



log 



+ -lo 



f-l0£ 



2 (Aa_i,Ac+l) 2 (Aa+2,Ac) 2 (Aa+2, Ac+2) 2 (Aa+2, Ac+1; 

The condition 



(Aa, Ac)(Aa+l5 Ac+2) 
(Aa+1, Ac)(Aa, Ac+2) 

must be true for rescaling invariance to hold. 
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(129) 
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